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Abstract 

Semiclassical expansions derived in the framework of the Extended Thomas- 
Fermi approach for the kinetic energy density r(r) and the spin-orbit density 
J(r) as functions of the local density p(r) are used to determine the central 
nuclear potentials V n {r) and V p {f) of the neutron and proton distribution for 
effective interactions of the Skyrme type. We demonstrate that the conver- 
gence of the resulting semiclassical expansions for these potentials is fast and 
that they reproduce quite accurately the corresponding Hartree-Fock average 
fields. 
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1 Introduction 

Mean-field calculations have been extremely successfull over the last 3 decades to 
describe the structure of stable as well as radioactive nuclei and this over a very 
wide range of nuclear masses. Especially effective nucleon-nucleon interactions of 
the type of Skyrme [|l], |] and Gogny forces |J have been particularly efficient in 
this context. Such phenomenological effective interactions can be understood as 
mathematically simple parametrisations of a density-dependent effective G-matrix 
(see H and || for a review on such effective forces). 

Together with the exact treatment of the mean-field problem in the Hartree-Fock 
(HF) approach, semiclassical approximations thereof have proven very appropriate. 
Especially the approach known as the Extended Thomas-Fermi (ETF) method has 
been shown to describe very accurately average nuclear properties in the sense 
of the Bethe-Weizsacker mass formula J7], In their selfconsistent version the 
ETF calculations determine the structure of a given nucleus by minimizing the 
total energy with respect to a variation of the neutron and proton densities. Such 
calculations require, however, only integrated quantities as e.g. the total nuclear 
energy. 

The aim of the present paper is, on the contrary, to investigate local quanti- 
ties such as nuclear mean-field potentials, effective mass and spin-orbit form factors 
which are at the basis of the description of the nuclear structure and which can 
be obtained as function of the selfconsistent semiclassical densities. The conver- 
gence of these local (non integrated) semiclassical quantities and their comparison 
to the corresponding HF distributions has, to our knowledge, never been extensively 
investigated as this will be done here. 

This paper is organized as follows. After specifying in section 2 the ETF ex- 
pressions up to order h 4 for r[p\ and J[p] using the general but rather cumbersome 
form of these expressions given in ||, [TO] , we show in section 3 that the semiclassical 
expansions which define these quantities converge rapidly for reasonable forms of 
the nuclear density p(r). Once these expressions and their convergence established, 
we compare in section 4 the average mean fields obtained using these semiclassical 
densities to the central potentials resulting from a HF calculation. We finally con- 
clude giving an outlook on how these calculations can be generalized to excited and 
rotating nuclei. 

2 Form factors for Skyrme Interactions 

For effective nucleon-nucleon interactions of the Skyrme type the total energy of a 
nucleus is a functional of the local densities p q (r), the kinetic energy densities r q (f) 
and the so called spin-orbit densities J q {r) 




(1) 
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where the subscript q = {n, p} denotes the nucleon charge state. In the case of 
broken time-reversal symmetry the energy density depends, in addition, on other 
local quantities |TI|, 0, such as the current density j(r) and the spin- vector density 



p(r). In what follows we will restrict ourselves to time- reversal symmetric systems 
leaving the case of broken time-reversal symmetry, particularly encountered in the 
case of rotating nuclei, to a subsequent publication. 

The total energy determined in this way is exact within the HF formalism. A 
semiclassical approximation thereof is obtained when instead of using the exact 
quantum-mechanical densities p q (r), T q (r) J q (r), etc. a semiclassical approximation 
for these quantities is used. The semiclassical densities r q {r) and J q {r) are obtained 
in the so called Extended Thomas-Fermi model as functions of the local density 
p(r) and of its derivatives. The best known of these functionals is the Thomas-Fermi 
approximation for the kinetic energy density in the form 

Tr ] {p«] = l^ 2 ) 2/3 pi /3 . <h^>- (2) 

Once these functional expression are given, the total energy of the nuclear system 
is then uniquely determined by the knowledge of the local densities of protons and 
neutrons. That such a functional dependence of the total energy on the local den- 
sities p q {r) should exist is guaranteed by the famous theorem by Hohenberg and 
Kohn fl3| . In the most general quantum mechanical case this functional is, how- 



ever, perfectly unknown. The great advantage of the semiclassical approach used 
here, consists in the fact that, in connection with effective interactions of the Skyrme 
type such an energy functional £ can be derived explicitly. In addition, it is to be 
noted that the semiclassical functionals obtained in the ETF formalism such as r[p] 
are completely general and valid for any local interaction and any nucleus, and can 
therefore be determined once and forever. 

Once the functional of the total energy is known, one is able, due to the Hohenberg- 
Kohn theorem, to perform density variational calculations, where the local densities 
p q (r) are the variational quantities. One should, however, keep in mind that, as the 
ETF functionals are of semiclassical nature, the density functions p q (r) obtained as 
a result of the variational procedure can only be semiclassical in nature, i.e. of the 
liquid-drop type. Taking into account that in such a process the particle numbers iV 
and Z should be conserved, one can formulate the variational principle in the form 

S p 

{£[p n (r),p p (r)] - X n p n {r) - X p p p {r)} d 3 r (3) 



with the Lagrange multipliers A n and X p to ensure the conservation of neutron and 
proton number. 

This density variational problem has been solved in two different ways in the 



past: either by resolving the Euler-Lagrange equation [15, [L6l resulting from eq. 



or by carrying out the variational calculation in an restricted subspace of functions 
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adapted to the problem, i.e. being of semiclassical nature, free of shell oscillations in 
the nuclear interior and falling off exponentially in the nuclear surface. It has been 



shown || [T5| that modified Fermi functions which for spherical symmetry take the 
form 

,r — Rq 



1 + exp(- 



-79 



q={n,p} 



(4) 



are particularly well suited in this context and that the semiclassical energies ob- 



tained are, indeed, very close to those resulting from the exact variation [15, 16 
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Fig. 1. Comparison of selfconsistent neutron and proton HF (solid line) and ETF 
densities (dashed line) for 208 Pb calculated with the SkM* Skyrme force. 

As an example of the quality of the semiclassical density obtained by such a 
restricted variation in the subspace of modified Fermi functions, eq. (01), we show 
in Fig. 1 a comparison of the neutron and proton densities obtained in this way 
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within the ETF approach with the corresponding Hartree-Fock (HF) densities for 
the nucleus 208 Pb calculated with the Skyrme interaction SkM* |TJ|]. It should be 
emphasized here that a similarly good agreement as the one shown in Fig. 1 is 
obtained for other nuclei or using other effective interactions such as the Skyrme 
forces SIII |7| and SLy4 |B . 

The energy density £ appearing in eqs. (0) and (fj) can be written for a Skyrme 



interaction as defined in ref. JT9[ in the compact form 



Sir) 



2m 



r + B lP 2 + B 2 (p 2 n + p 2 p ) + B 3 pr + B 4 (p n r n + p p r p ) 



B 5 (Vp) 2 -B 6 {V Pn f + {V Pp f + p a [B 7 p 2 + B 8 (p 2 n + p 2 )} 



-Bo 



J-Vp + J n - Vp n + J p ■ Vp p + S Co ui(r) , 



(5) 



given in terms of the coefficients Bi — Bq (see table 1) instead of the usual Skyrme 
force parameters t , ti,t 2 , t 3 , x , xi, x 2 , x 3 , W : 
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Tab. 1. Correspondence between the coefficients Bi used in the text and the usual 
Skyrme force parameters. 



In eq. (|5|) non-indexed quantities like p correspond to the sum of neutron and proton 
densities as p— p n + p v and Seoul is the Coulomb energy density which can be written 
as the sum of the direct and an exchange contribution, the latter being taken into 



account in the Slater approximation J2(| |21 



Scouiif) = ^ppir) I d 3 r' 



.Pp^X _ 3 2(3)1/3 p 4/3 M 
I r — r 1 1 4 7T 



(6) 



The HF equation is obtained through the variational principle which states that 
the total energy of eq. (P should be stationary with respect to any variation of the 
single-particle wave functions : 



-v 



h 2 

2m* (r) 



V + V q (r) - iW q {r) ■ (V x a 




(7) 
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Here appear different form factors such as the central one-body potential V q (f), the 
effective mass m*(r) and the spin-orbit potential W q (r) which are all defined as 
functional derivatives of the total energy density. One obtains from eq. (H) : 



SSjf) 

2(B lP + B 2 p q ) + B 3 t + B 4 r q + 2(B 5 Ap + B 6 Ap„) + (2 + a)B 7 p 



a+l 



+B 8 



+ i?9(divJ + divJ 9 ) + Vcouiir) 5. 



pq ) 



(8) 



and 



m 



W q {r) 



2m 5S(f) 
~W 6r q (f) 

5£(r) 



Off) 

1 + ^ [B 3P (r) + B, Pq (f)] 



-Bg V{p + Pq ) . 



The Coulomb potential in eq. (§) is easily obtained as 



(9) 



(10) 



fill 



7T 



and the spin- 



It is noteworthy in this connection that for such an effective mass 
orbit potential flTOP the energy density (^) takes the simple form 

= ^Y,fq^q + B lP 2 + B 2 {pl + pi) - B 5 (Vp) 2 - B 6 [(Vp n ) 2 + (Vp 



+p a [B 7 p 2 + B 8 (pl + pi)} + £ J q ■ W q + Scouii?) (12) 



which simplifies somewhat the calculation. 

All the expressions derived so far (eqs. (||) - ([12])) are exact and when used with 
densities constructed from the single-particle wave functions ip ( J 9 \f), solutions of the 
HF equation ([?]) these quantities contain all the quantum effects of the system. If 
one is interested in the semiclassical approximation of these form factors one can 
immediately conclude, from the analytical form of eqs. (|9]) and (|1(]) and the smooth 
behavior of the semiclassical densities, as demonstrated in Fig. 1, on the smooth 
behavior of the semiclassical effective mass and spin-orbit form factors. As the 
nuclear quantal density is well reproduced on the average by the ETF densities it 
appears evident that the same is going to be the case for the effective mass and the 
spin-orbit potential, when semiclassical, i.e. liquid-drop type densities are used in 
eqs. © and 
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Things are, however, less evident for the central nuclear potentials. Not only is 
the functional expression, eq. (|8|), much more complicated then those for the effective 
mass and spin-orbit potential, but the central potential is also the only of the three 
functional derivatives that depends not only on the local densities p q (r) and their 
derivatives but also on the kinetic energy density T q (f) and the spin-orbit density 
J q {r) which are the quantities for which the Extended Thomas- Fermi approach has 
written down functional expressions. We therefore choose to study the convergence 
of the semiclassical series corresponding to these functional expressions of r q [p q ] 
and Jq[p q ] before investigating the quality of the agreement between the HF central 
potential and the one obtained when using these semiclassical ETF functionals. 



3 Convergence of ETF functionals 

The semiclassical expansions of kinetic energy density r q and spin-orbit density J q 
as functions of the local density p q are functional expressions with h as order pa- 
rameter. These expressions can be obtained for instance through the semiclassical h 
expansions developed by Wigner El and Kirkwood EM or through the semiclassical 



method of Baraff and Borowitz p5| . In either of the two approaches one obtains 
functional expressions like 

i ETF) = \p*] + \p*\ + ^ \p,] (13) 



written here for the kinetic energy density r q {r) where T^ TF '[p q ] is the well known 
Thomas- Fermi expression already given in eq. (0), T^[p q ] the semiclassical correc- 
tion of second order and [p q ] is of fourth order in h. The ETF expressions such 
as rf TF ) [p q ], eq. (|l3l) , up to order h are to be understood as the converging part 
of an asymptotic series. 

The second order term [p q ] has already been derived in ref. Q for a Hamil- 



ivp mass rr u ^ 

2 



tonian, eq. (|7|), with an effective mass m* = mj f q and a spin-orbit potential W q 



T. 



1 (V Pg ) 2 1 A 1 Vpq ■ Vf q , 1 Af q 1 Vf q 



36 Pq 3 ry 6 f q 6 rq f q 12 ry \ /, 

1 (2m\ 2 (W q \ 
+ 2 iw) "« (tj (14 > 

the first term of which is known as the Weizsacker correction 0. It was sometimes 
used in the past with an adjustable parameter (instead of ^) in order to mock up 
the absence of the other second order and all of the fourth order terms. It has, 
however, been shown (see e.g. ref. ^) that such a procedure is unable to correctly 
describe both the slope of the surface of the nuclear mass or charge density and, at 



7 



the same time, the height of nuclear fission barriers in the actinide region. From 
this analysis we conclude that the inclusion of fourth order terms in the semiclassical 
functionals is, in fact, without credible alternatives. 

In the following we are going to exploit the expressions for the 4 th order function- 
als T q and J q developed by Grammaticos and Voros ||, [K| . These authors have taken 
the convention "that any free-standing gradient operator acts only on the rightmost 
term" (see their remark after eq. (III. 13) of ref. |J). In our present work we prefer 
to rewrite these terms in a more conventional way and have any free standing gra- 
dient operator act, as usual, on all the terms that appear on its right hand side. We 
therefore write (subscripts GV refer to the Grammaticos/ Voros convention) 



(W 9 -v) 2 /, 



GV 



;(V/,-V)(V/, 



(Vp 9 -V) 2 / 9 



GV 



Vp,-V(V/,-Vp,) - -VL-V{V Pq f 



and 



{Vf q -V){Vp q -V)L 



GV 



;V/yV(V/ 3 



plus terms that are obtained from these ones by interchanging the role of f q and p q . 
One has also to keep in mind that Grammaticos and Voros use a slightly different 
definition of the effective-mass and spin-orbit form factors than the ones given in 
eqs. @ and (|10D above : 

J GV J 



m 



and 



s. 



Using these expressions we obtain the following form for the 4 th order kinetic 
energy density, where contributions from the spin-orbit interaction have been, tem- 
porarily left out. 
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fqPq 



fqPq 
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- 16 
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f 2 P 2 

Jq Hq 



fqfi 



(15) 



The interested reader, who might want to use the semiclassical functionals calcu- 
lated here, will find in the Appendix the expression that the 4 th order kinetic energy 
density takes in the case of spherical symmetry as well as all the other semiclassical 
functionals developed below. 

Until now we have not taken into account the spin-orbit interaction. Its influence 
on the semiclassical ETF functionals is treated in ref. |10| and its contribution r^ so 
constitutes simply an additive term to the spin-orbit independent part of the kinetic 
energy density considered above. According to JTOl] : 



T, 



(4) s 



2 n l/3 



[p] = (3. 2 )- 2 /3 (%) ± 



UV 4/2 



-w q -Aw q + -wvv(divwg + -a(w c 



+ -(dwW g 



2W q ■ {Vf q -V)W q + divW, {Vf q -W q ) + Vf q ■ (W q - V)W q 



+W 2 q Af q + W q -V(W q -Vf q ) - -Vf q -V(Wl 



1 



+ 



4/ 2 



^fqfW 2 q +{W q -Vf q 



Qpq 



W q -{V Pq -V)W q + (W q -Vp q ) divW q 



SfqPq 



{Vf q -V Pq )W 2 q + {Vf q -W q ){V Pq -W q ) 



(16) 



It is now interesting to investigate the relative importance of the different contri- 
butions in eqs. (@), fll4T) , (p~5|) and dl6|) to the kinetic energy density obtained when 
using the selfconsistent semiclassical densities generated by a variational procedure 
restricted to functions of the type of eq. (|j) as explained above. As can be seen on 
Fig. 2 (a) the Thomas- Fermi contribution to r[p] is largely dominant, at least in the 
nuclear bulk. Semiclassical corrections play, however, a significant role in the nuclear 
surface with a second-order correction which is much larger than the fourth-order 
term (multiplied for better visibility by a factor 10 in Fig. 2 (a)). The different con- 
tributions to the second and fourth order functional are given respectively in part 
(b) and (c) of the figure. We show the contributions coming form gradient terms of 
p (term 1 and 2 in eq. (|T4T)), of / (terms 4 and 5) and the mixed term (term 3) as 
well as the spin-orbit contribution (last term) and similarly in part (c) of the figure 
for the fourth-order term. As can be seen, the gradient term of p is dominant in 2 nd 
order, whereas in 4 nd order the spin-orbit contribution becomes also crucial. The 
selfconsistent HF neutron kinetic energy density is also shown on Fig. 2 (a). One 
notices that, except for quantum oscillations in the nuclear interior, the HF kinetic 
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energy density is quite nicely reproduced if the semiclassical corrections T2 and T4 
are taken into account. 
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Fig. 2 Contributions from the different orders in the semiclassical expansion to the 
kinetic-energy density r[p] for the self consistent neutron density distribution shown 
in Fig. 1 for 208 Pb (TF (solid line), 2 nd (dashed line) and 4 th -order multiplied 
(dotted line)) are compared with the corresponding HF density (dash-dotted line) 
(part (a)). Different contributions to 2 nd (part (b)) and 4 th order (part (c)) coming 
from gradient terms of p (solid line), of / (dashed), of mixed terms containing 
gradient terms of p and / (dotted), and of the spin-orbit coupling (dash-dotted 
line) . 
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It is interesting to note in this connection that despite the fact that the 2 nd -order 
contribution (f) is one order of magnitude larger than the 4 th -order contribution, 
after integration the 2 nd -order contribution of r[p] to the total energy, i.e. the inte- 
gral J2 q I fq(r)Tq 2 K'^)d 3 r is of the same order of magnitude than the corresponding 
4 th -order contribution (see e.g. [57, M), which seems to indicate that there is a 
stronger cancellation taking place in the 2 nd -order than in 4 th -order contribution. 

We have done the same study for the proton distribution, for other nuclei and 
used other effective interactions of the Skyrme type, namely the Skyrme SIII force 
H[] and the SLy4 Skyrme force ||18|| . The conclusions made above remain valid in 
all these cases, only the relative importance of the semiclassical corrections [p] 
and r^fp] increases slightly when one goes from heavy to light nuclei. 

As already discussed above, ETF functionals like r^ ETF ^[p], eq. (0), up to or- 
der h 4 constitute the converging part of an asymptotic expansion which needs to 
be truncated. Comparing, indeed, the p dependence of the different orders of the 
semiclassical functional (p 5//3 for r^ TF '[p\, p for r^fp] and p 1 / 3 for r^[p]) one con- 
cludes that a term [p] in the ETF functional would show a p dependence of the 
form p -1 / 3 and would therefore diverge in the limit r — > oo for densities that fall off 
exponentially at large distances. 



Let us now turn to the spin-orbit density J. It is given in ref. [|TU[ in the form 
of second rank tensor which is related to the components of the vector J by the 
relation 

flU 

where e\^ v is the Levi-Civita symbol. The spin being a purely quantal property with 
no classical analogon, there is no contribution to the semiclassical functional of J in 
lowest order, i.e. at the level of the Thomas-Fermi approach whereas one obtains for 
the 2 nd and 4 th -order contributions to the semiclassical expansion of the spin-orbit 
density 

/(2) 2mp, # 
h 2 f q Wq 

and 



(17) 
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(V Pq -V)W q + divW q V Pq 



f 2 

Jq 
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'7a 



W q Af q + (W q -V)Vf q 

2m, 



W q ) 2 w q + (vf q -w q )vf q - 2 (=^y w q 



Vf q -VPa)W q +(Vf q -W q )V Pc 



Let us again investigate the convergence of the semiclassical expansion associated 
this time with the vector field J q [p] and compare it with the corresponding HF 
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density. We show in Fig. 3 the contributions to J q [p] from 2 nd and from 4 th order 
as well as Jhf- 



0.01 — 







Pb 






\ SkM* 








I I T I T I T I 1 



r (fm) 

Fig. 3. Contributions from the different orders in the semiclassical expansion (2 nd 
(solid line) and 4 th -order (dashed line)) to the radial part of the vector field J q (r) 
shown here for the selfconsistent spherical neutron distribution of 208 Pb as compared 
to the HF spin-orbit density (dash-dotted line). 

We would like to check now that the semiclassical functionals which we have 
written down up to order h 4 are indeed correct. We perform this test numerically 
in the following way : 

One notices that when calculating the total energy through eq. ( |T2"D the kinetic 
energy density r q does not appear by itself, but only in connection with the form 
factor f q . The r-dependent part of the total energy is simply obtained through the 
integral / d 3 rJ2 q fq T q- I n the contribution at order h 4 to this integral one can then 
perform integrations by parts to obtain an expression which contains only second 
order derivatives of the density p q and the effective-mass form factor f q . 

d 3 r- 




3240 y p q ) Pq 
where t^ 4 ' is the spin-independent part of the kinetic energy density. 
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The dependence of the total energy density on the spin degrees of freedom en- 
ters in two different ways : through the spin-orbit part of the kinetic energy den- 
sity Tg° and through a term J q -W q . This total spin-dependence is then given by 

^9 J (^>mf<i T q [p\ ^9 ' ( ^ r - ^ can be sriown after some integration by parts 

with eqs. fllij) , ([TBI), dO) an< ^ dUD that i* 1 the different orders of the semiclassical 
expansion this integral takes on the simple form B 



and that 



— fA A Hp\ + w q -j^[ P ] 




(20) 



(21) 



We have tested the semiclassical functionals given above by verifying numerically 
that these integral relations (|19|) , (|20| ) and ( pl|) hold true. 

We have also evaluated different integrals involving these functionals and which 
have been calculated in ref. (jlfj. We obtain agreement with their results of the 
order of 1 to 2 %, which is of the same order as their agreement between the results 
of a full variational calculation and one in the restricted subspace of modified Fermi 
functions. 

As can be seen on eq. (|||) only the divergence of the vector field J is present in 
the expression of the central one-body potential. One obtains from eqs. ( |T7D and 
respectively the contributions to 2 nd order 



Pq dwW q + (V Pq -W q ) - ^(Vf q -W q ) 

Jq 



[22) 



and, after some lengthy but straightforward calculation, to 4 order 
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+3A/,(V/ 9 -W,) + 5Vf q -(W q -V)Vf q + bVf q -(Vf q -V)W q 
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(23) 



4 Average nuclear potentials 



One now controls all the ingredients which enter into the calculation of the nuclear 
central potentials, eq. (d), for effective interactions of the Skyrme type. It is now 
interesting to look at the convergence of the expressions which define these average 
fields calculated with the semiclassical functionals r q [p] and J q [p] and to check how 
these potentials compare with the ones obtained in the HF approach. 

We therefore show in Fig. 4 the neutron and proton nuclear central potentials 
for the nucleus 208 Pb obtained for the Skyrme interaction SkM*. The Coulomb 
potential Vc ou i for the proton field has been left out in this investigation, because it 
is directly given through the proton density (see eq. flnD). As the latter is very well 
reproduced, except for quantum oscillations in the nuclear interior, we already know 
that the Coulomb potential calculated through this semiclassical density will, indeed, 
reproduce on the average the exact one calculated from the quantum-mechanical 
densities. 

We show a comparison between the HF neutron and proton central potentials 
with the ones obtained using the selfconsistent semiclassical densities p n and Pp but 
restricting ourselves to the TF approximations for the functional r[p], eq. (pj), and 
J[p] (which is zero as explained above). We do not want to call this the Thomas- 
Fermi approximation to the nuclear central fields since even if we have used the above 
mentioned functionals in their Thomas-Fermi approximation, the nuclear structure 
has been determined through a full variational calculation including the functionals 
up to order h 4 . 

As seen on the figure the reproduction of the HF selfconsistent fields is already 
quite remarkable at the lowest (TF) order in the semiclassical expansion. Apart 
from shell oscillations in the nuclear interior and a small wiggle in the TF potential 
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in the surface region the agreement seems very satisfactory. 



o 



SkM" 



r (f m ) 




O 2 4 6 8 10 12 

r (fm) 



Fig. 4. Comparison of the Hartree-Fock central nuclear potentials (solid line) for 
protons and neutrons with the corresponding semiclassical potentials (dashed line) 
obtained using the Thomas- Fermi approximation for the functionals r[p] and J[p\. 

It is now interesting to study the contributions to the nuclear central fields 
coming from higher orders in the semiclassical expansion. For this reason we show 
in Fig. 5 the corrections SVjf 1 and 5V^ denned as (see eq. (§)) 

6V™ = (B 3 + BJt® + B*r® + B 9 (2 divjf + divjf ) , 

5K (4) = (Bs + B 4 )t^ + S s r« + B 9 (2 div JJ 4 > + divjj 4 >) (24) 

together with the semiclassical TF potential already shown in Fig. 4. It turns out 
that these corrections are rather small and we have to multiply 5V^ by a factor of 
10 and 8V^ by a factor of 100 to make their relative importance better visible on 
figure 5. 
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Fig. 5. Semiclassical neutron potential using the TF approximation (solid line) 
of the semiclassical functionals r[p], and J[p] and corrections coming from second 
(dashed) and fourth order (dotted line) in the semiclassical expansion. For better 
visibility the second order correction has been multiplied by a factor of 10 and the 
fourth order by a factor of 100. 

It can be seen that both these terms give a contribution in the nuclear surface 
where the lowest order term showed some deviation from the HF potentials. It is 
therefore to be expected that potential using the second order functionals [p] and 
[p] will partially correct for this deficiency and be quite close to the selfconsistent 
HF potentials. Due to the smallness of the 4 th order term in Fig. 5 we can expect 
the semiclassical potentials obtained using the full semiclassical functionals up to 4 th 
order to be practically indistinguishable from the ones using the 2 nd order corrections 
only. This conclusion is, indeed, confirmed on Fig. 6. 
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Fig. 6. Comparison of the Hartree-Fock central nuclear potentials (solid line) for 
protons and neutrons with the semiclassical potentials (dashed line) obtained by 
using the semiclassical functionals r[p] and J[p] up to order h 2 (dashed line). The 
semiclassical potentials including the functionals up to order h 4 are indistinguishable 
from the latter ones. 

The same kind of calculations have been performed also for lighter nuclei down 
to 40 Ca. The quality of the agreement between HF and semiclassical potentials 
is the same as the one obtained for the nucleus 208 Pb studied above. As already 
mentioned, the effective-mass form factor f q (r), eq. (|9[) and the spin-orbit potential 
W q (r), eq. ATI]) are, except for shell oscillations in the nuclear interior, very well 
reproduced, since the nuclear densities which directly determine these quantities are 
well reproduced. For this reason we do not explicitly show these quantities. These 
conclusions remain valid when other effective interactions of the Skyrme type are 
used. 

5 Summary and conclusions 

We have demonstrated that using the Extended Thomas-Fermi approach one is not 
only able to give a very precise description of average nuclear properties but that 
this method is also able to reproduce quite nicely local quantities, not only neu- 
tron and proton density distributions but also the corresponding nuclear central 
potentials, effective-mass form factors and spin-orbit potentials. These are precisely 
the ingredients of the Schrodinger like Hartree-Fock equation, eq. fl7|), which arises 
from the variational principle. Within this semiclassical approach which relies on a 
density variational calculation one should therefore be able to solve in an approxi- 
mate way the quantum mechanical problem without having to go through the full 
selfconsistency problem of the HF approach. This is the essential idea of an ap- 
proximate solution of the HF problem known as the "expectation value method" 
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(EVM) p8| . It consists in constructing the ground state Slater determinant from 
the eigenf unctions of eq. (0) using the ETF fields V q (f), m*(r) and W q {r) (to second 
order in the ETF expansions) and calculating with this Slater determinant the ex- 
pectation value of the total Skyrme Hamiltonian. The rapid convergence of the ETF 
functionals demonstrated above explains, a posteriori, the success of this approach. 

The interested reader might object that nowadays, where computational power 
has been increased tremendously, there is no real need for semiclassical approxima- 
tions, but all calculations of nuclear structure should be directly performed at the 
level of the Hartree-Fock model (or beyond). The point, however, is that as soon as 
one is interested in nuclear systematics where one is looking at the behavior of nuclei 
over a wide range of the nuclear chart, semiclassical approximations are without any 
credible substitute. 

A point that might e.g. be interesting to study is the variation of the diffuseness 
of the nuclear densities and central potentials when increasing the nuclear excitation 
and/or when going to rotating nuclei. The approach we have developed here is, 
indeed, easily generalized to the description of excited or rotating nuclei. If one 
is interested in hot nuclear systems one simply need to replace the semiclassical 



functionals for r[p] and J[p] which we have developed by those derived in ref. [R 
for the ETF approach at finite temperature. In this case the coefficients in the 
semiclassical expansions, eqs. (0), (|HD, ([[§), ( |ToD etc. are to be replaced by a 
combination of Fermi integrals |15| . 

If on the other hand one is interested in the description of systems breaking time 



reversal symmetry, the energy density, eq. (|5|), will be changed |TT], |IIJ and some 
additional densities will appear, such as the current density j(r) or the spin- vector 
density p(r) which, in the case of rotations, are a manifestation of the time-odd 
part of the density matrix generated by the cranking piece of the Hamiltonian. This 
causes not only a change in the analytical form of quantities such as the average po- 
tentials V q {f) which are now going to depend on these additional densities but also 
leads to the appearance of additional terms in the functionals r[p] and J[p\. These 
functionals have already been determined in \Y1\ which makes it quite straightfor- 
ward to calculate selfconsistent semiclassical fields such as V q (r) for rotating nuclei 
very similarly to what we have done here in the static case and investigate the 
dependence of these quantities with increasing angular momentum. Investigations 
along these lines are currently in progress. 

The present method can also be profitably exploited to establish some easy to use 
parametrization (as function of mass number A and isospin parameter I=(N-Z)/A) 
of central and spin-orbit potential and effective mass and this over a wide region of 



the nuclear chart ||29|| . Such an investigation could be the ideal starting point for 
Strutinski shell correction calculations as e.g. formulated in the so-called Extended 
Thomas- Fermi plus Strutinski integral (ETFSI) method (see ref. [3(J and references 
given therein) where average mean fields like the ones investigated here are used to 
determine shell corrections. 
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Appendix 

As mentioned in the text we gather here the expressions the ETF functionals 
take in the case of spherical symmetry. 

For the 4 th order spin independent part of the kinetic energy density, eq. (|15D , 
one obtains (primes denoting derivatives with respect to the radial variable r) : 
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When giving the spin dependent part of T^ 4 )[p] we take advantage of the fact 

that the spin-orbit potential W q has for Skyrme forces the simple form of eq. flTD| ) 
which allows us to introduce the quantity 

A g = p + Pq W q = -B 9 VA q (A.2) 

Using the vector identity 

rot (rot a) = V(div a) — A a 

one shows that because of the form of the spin-orbit potential, eq. (|i~0D , one simply 
has 

W-AW = W-V(divW) 

which then allows us to write the spin dependent part of r^ s ° in the form of the 
local densities p n and p p and of their derivatives : 
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which for spherical symmetry reads 
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The 4 th order spin-orbit density J^\p\, eq. (18) written in terms of the function 



A q defined above is given by 
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which in the case of a spherically symmetric system takes the form 
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with the unit vector in radial direction u r . 

The corresponding expressions for div [p] are the following : 
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and in the case of spherical symmetry 
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which can also be obtained directly from eq. (|A.6 ), remembering that the divergence 
of a vector field a that has only a radial component a r is given by 

diva = \-^-(r 2 a r ) . 
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